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Abstract 

We study the dynamical formation of disoriented chiral condensates in very high energy 
nucleus-nucleus collisions using Bjorken hydrodynamics and relativistic nucleation theory. 
It is the dynamics of the first order confinement phase transition which controls the 
evolution of the system. Every bubble or fluctuation of the new, hadronic, phase obtains 
its own chiral condensate with a probability determined by the Boltzmann weight of the 
finite temperature effective potential of the linear sigma model. We evaluate domain size 
and chiral angle distributions, which can be used as initial conditions for the solution of 
semiclassical field equations. 



1 Introduction 



It is expected that quark-gluon plasma will be created in high energy heavy ion collisions 
at the Relativisitic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory and 
at the Large Hadron Collider (LHC) at CERN. The experiments should allow us to study 
the chiral and/or confinement phase transition/crossover. One possible consequence of 
plasma formation and disassembly is the creation of misaligned chiral domains (disoriented 
chiral condensates or DCC) which are regions of space-time where the average value of 
the chiral field does not point in the same direction as the surrounding true vaccuum. For 
a brief review see [|1|. 

This possibility was really emphasized and studied by Rajagopal and Wilczek || in 
the context of heavy ion collisions within a version of the linear a-model. It was found 
that if the system evolves close to thermal equilibrium the average domain size is of order 
1/T C , the inverse of the phase transition temperature. The resulting domain sizes are 
therefore too small to have many observable consequences; enhanced baryon-antibaryon 
production may, however, be one of them ||. Subsequent studies of DCC concentrated 
on nonequilibrium scenarios like quenching 0] where the system is suddenly relaxed from 
an initial thermal state above T c to zero temperature. The resulting field configuration is 
unstable and decays on time scales of order l/m a . Again one obtains domain sizes which 
are too small to have many experimental consequences. Finally, in an annealing scenario 
||, medium modifications of the a were taken into account, allowing for a reduction of 
m a and therefore an increase of domain size of up to 3 fm. 

All aforementioned approaches assume that the phase transition is, or at least close 
to, second order. In this paper we would like to investigate the consequences of a first 
order phase transition. The disoriented chiral domains are created via statistical fluctu- 
ations. We assume that a relativistic heavy ion collision results in the formation of an 
extended volume of quark gluon plasma which transforms itself into a hadronic resonance 
gas through homogeneous nucleation of hadronic bubbles Bubbles of the new 

hadronic phase are created by statistical fluctuations which grow due to their lower free 
energy. The phase transition is completed when the entire plasma has been transformed 
into the hadronic resonance gas. This dynamics implements the confinement character- 
istics of the phase transition. The dynamical information is contained in the free energy 
difference between the plasma phase and the resonance gas phase and in their transport 
properties. 

We assume that the dynamical evolution will be dominated by this color confinement 
at temperatures close to the critical temperature T c . In principle we could imagine that 
every hadronic bubble nucleated in this way starts out with the chiral field pointing in 
a random direction in internal symmetry space. The bubbles grow with time until the 
whole plasma phase is transformed into the hadronic phase, resulting in a distribution of 
domains with different sizes and different chiral angles. For definiteness we use the linear 
cr-model to evalutate the part of the hadronic free energy dependent on the chiral field. 



2 



The effect of this contribution on the dynamical evolution is perturbative. 

The domains formed in this approach will be larger or of order R c , the critical size 
radius to nucleate an hadronic bubble. This is a dynamical quantity which diverges at the 
critical temperature and approaches the hadronic length scale 1/T C at small temperatures. 
We calculate a domain size and chiral angle distribution. This distribution may be taken 
as the initial condition for solving the semiclassical equations of motion for the subsequent 
time evolution of the system ||. 

Our approach will neglect the possibility that a bubble nucleated with a very large 
radius might contain several chiral domains. This is not a serious concern since the 
average nucleation bubble is not so large. Furthermore, the resulting domain size and 
angle distribution should also depend on the rapidity. For given rapidity the nucleated 
bubbles will obtain a certain eccentricity in the beam direction. In this first study we 
neglect this complication. 

In section 2 we start out by outlining the implementation of the dynamical evolution 
of the system. The free energies we are using in this section are at zero chiral angle. In 
section 3 we evaluate the perturbative chiral contribution to the hadronic free energy at 
nonzero chiral angle and show how to obtain the domain size and angle distribution. In 
section 5 we discuss and summarize our results. 

2 Dynamics of the Confinement Phase Transition 

In this section we briefly review the dynamics of the hot matter as it passes through a 
first order confinement phase transition. The overall picture is that quark-gluon plasma 
is formed at high temperature and subsequently expands and cools. The plasma must 
cool below the critical temperature before bubbles of the hadronic phase can be nucleated. 
When enough bubbles have been nucleated their growth causes a reheating, turning off 
further bubble nucleation. Eventually all matter is converted to the hadronic form, and 
cooling begins again. Eventually the hadrons lose thermal contact with each other and 
they free-stream to the detectors. It is within this scenario that we consider the formation 
of DCC. The works should be consulted for more details about this scenario. 

2.1 Review of nucleation dynamics 

The rate for the nucleation of the hadron phase out of the plasma phase can be written 

as 

I = I e- AF */ T , (1) 

where AF* is the change in the free energy of the system with the formation of a critical 
size hadronic bubble and Jo is the prefactor with dimensions of inverse volume inverse 
time. In general, statistical fluctuations at T < T c will produce bubbles with associated 
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free energy 

Air 

AF = T [p q (T) - p h (T)]R 3 + AnR 2 a. (2) 

Here p is the pressure of the quark or hadron phase at temperature T, and a is the surface 
free energy of the quark- gluon/hadron interface. Since p q < p^ it follows that there is a 
bubble of critical radius 

R4T) = 2a/\p h (T)-p q (T)}. (3) 

Smaller bubbles tend to shrink because the surface energy is too large relative to volume 
energy, and larger bubbles tend to grow. The free energy of the critical size bubble is 
therefore 

AF* = ^kuRI (4) 

The prefactor has been computed in a coarse-grained effective field theory approximation 
to QCD to be 

3tt \3TJ ^(Aw) 2 ' 1 ' 

where r] q is the shear viscosity in the plasma phase, £ q is a correlation length in the plasma 
phase, and Aw is the difference in the enthalpy densities of the two phases. At the critical 
temperature, R* — > oo, and the rate vanishes. 

Given the nucleation rate one would like to know the (volume) fraction of space h(t) 
which has been converted from QCD plasma to hadronic gas at the proper time t, which is 
the time as measured in the local comoving frame of an expanding system. This requires 
a kinetic equation which uses / as an input. If the system cools to T c at time t c then at 
some later time t the fraction of space which has been converted to hadronic gas is 

h(t) = f dt'I(T(t'))[l - h(t')]V buh (t',t). (6) 
Jt c 

Here Vb u b(t', t) is the volume of a bubble at time t which had been nucleated at the earlier 
time t'. 

Once formed, bubbles will grow, as it is favorable from the point of view of free energy. 
We assume the growth law 

v(T)=v [l-T/T c f 2 , (7) 

where t>o is a model-dependent constant. The simple illustrative model for bubble growth 
is 

V huh (t',t) = ^ (i*.(T(f)) + £ dt"v(T(t"))^ . (8) 

This expression assumes that the interface between the inside and outside of the bubble 
is created at rest in the local comoving frame. 
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A dynamical equation is needed to describe how the system expands. We will use 
Bjorken's longitudinal scaling hydrodynamics [|IIJ]. The derivative of the energy density 
is related to the enthalpy density as 

- - -- (9) 
dt t ' 1 } 

This assumes kinetic equilibrium among the particles but not phase equilibrium. It is a 
statement of energy conservation. The energy density is 

e(T) = h(t)e h (T) + [1 - h(t)]e q (T) 7 (10) 

where e^iT) and e q (T) are the energy densities in the two phases at temperature T, and 
similarly for w. 

We model the plasma phase by a gas of gluons and massless quarks of either two or 
three flavors with a bag constant B to simulate confinement. The pressure is 

p q = g^T*-B, (ii) 

where g q is chosen appropriately. We model the pressure in the hadronic phase according 
to the formula 

Ph = g^T (12) 

where is an effective number of degrees of freedom relevant for the temperature range of 
interest, namely, 130 < T < 170 MeV. Including noninteracting pions alone, for example, 
gives #h ~ 3. Here we will consider two parameter sets. Set A includes only u and d 
quarks. In the hadron phase we include the particles 7T, 77, 00, p. In the aforementioned 
temperature range ~ 4.6 |TTJ. Set B includes u, d and s quarks. In the hadron phase we 



include all states in the Particle Data Book with a repulsive interaction among them. This 



equation of state was studied in ||12|| . A good parametrization is obtained with « 7.5. 
The critical temperature is obtained by equating the pressures of the two phase. We 
adjust the bag constant so that T c = 160 MeV. In set A one gets B lj/4: = 220 MeV and in 
set B one gets 232 MeV. The other parameters are as in ||, namely: a = 50 MeV/fm 2 , 
£ q = 0.7 fm, rj q = 14.4T 3 and v = 3. 



2.2 Numerical results 

The evolution equations may be solved as in |J. Figure 1 shows the temperature as a 
function of time. It is assumed that the system cools to T c at a local comoving time of 
3 fm/c. The system must supercool to about 0.95T C before noticeable nucleation begins. 
The system continues to expand and cool due to inertia, but is slowed down by bubble 
nucleation and growth. The temperature reaches a local minimum at around 0.80T C , at 
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which point the release of latent heat is sufficient to begin reheating. As the temperature 
goes up, the nucleation rate decreases, and basically shuts off at about 0.95T C . Thereafter 
there is essentially no new bubble creation; the transition completes due to the growth 
of existing bubbles. When the fraction of space h occupied by hadrons reaches 100% the 
transition is complete and the temperature will again fall. This is not shown on the figure. 
It is the state of the system at this moment which will determine the initial conditions 
for solving the equations of motion for the chiral fields. 

Figure 2 shows the volume of a hadronic bubble at the end of the phase transition 
as a function of the time at which it was nucleated. The volume of a critical size bubble 
is a monotonically decreasing function of temperature. This means that bubbles which 
were nucleated early or late in the supercooling phase should have the largest volume, 
as is seen to be the case. The smallest bubbles should be those which were nucleated 
near the minimum of the the temperature curve, which is also seen to be the case. The 
minimum of the volume curve is offset to somewhat later times than the minimum of the 
temperature curve due to the subsequent growth of bubbles. 

Figure 3 shows the density of bubbles as a function of their volume at the end of 
the transition. This is obtained by keeping track of how many bubbles were nucleated 
between times t and t + dt and then following their growth to the end of the transition 
at time tf. The density distribution is 

^ = J t t ^I(t)[l-h(t)]5(V-V huh (t,t i )) , (13) 

where the factor t/tf takes into account the fact that the density decreases due to the 
expansion of the system. In general there will be two times contributing to the same 
final bubble volume, as already seen in figure 2. The distribution plotted in figure 3 
is made dimensionless by multiplying by the minimum volume squared and is plotted 
against the linear size of the bubble to better display the results. For set A the minimum 
final volume is 22 fm 3 while for set B it is 10 fm 3 . It should be remarked that at the 
end of the transition the bubbles will not be spherical due to the longitudinal, Bjorken, 
expansion along the beam axis. They will be elongated along the beam direction with a 
typical eccentricity of tf/t sc ~ 3 — 4 where t sc is the time at which the temperature has 
dropped to its minimum value along the supercooling curve. 

Some important numbers from the supercooling and reheating cycle are listed in Table 
1. For each paramter set we list: the time t sc and temperature T sc when supercooling ends 
and reheating begins; the completion time tf and temperature Tf of the phase transition; 
the minimum volume of a bubble at the end of the transition and the time t min at 
which it was nucleated. 

The distribution dn/dV at the end of the transition diverges at its lower limit like 

1 / \A ~ Vmin/V • This is a rather weak, integrable, singularity which may, however, have 
important consequences for the observation of DCC in heavy ion experiments. 

One remarkable feature of this dynamics is that the contribution to large bubbles 
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comes mainly from early times, t < t sc . Bubbles which are nucleated near the bottom 
of the temperature curve do not have sufficient time to grow very much before reheating 
causes the growth velocity to decrease substantially. Bubbles nucleated at later times 
have a larger (critical) volume to begin with, but their growth is also suppressed by a 
decreasing growth velocity. In addition, their rate of production is suppressed relative to 
early times due to the factor 1 — h in eq. (13); there is insufficient space available between 
existing bubbles to nucleate new ones. The distribution for volumes larger than about 
100 fm 3 is almost completely determined by the initial stages of the phase transition. 
Even if we would shut-off nucleation by hand at time t sc we would essentially obtain the 
same size distribution at large volumes. The reheating phase of the transition does not 
influence the size distribution at large volumes! 

The behavior of dn/dV at large volumes is well fitted to the form exp[— a (V/V min ) 2 / 3 ]. 
The values of the coefficient a are given in Table 1. The exponential dependence on the 
surface area is straightforward. Large bubbles at the end of the transition arise from 
large bubbles nucleated early on. The probability to nucleate one is proportional to the 
rate obtained from eqs. (1) and (4), which is exponential in the surface area of a critical 
bubble. Allowing for growth due to longitudinal expansion we can estimate roughly that 

2/3 2 

Knin- (14) 

If we compare this result with the fitted values given in Table 1 we get agreement to 
better than a factor of 2. 





/ 3 a 
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3 DCC Within a Nucleated Bubble 

In the previous section we discussed how a first order phase transition from the deconfined 
quark-gluon phase to the confined hadronic phase may occur via supercooling, nucleation, 
bubble growth, and reheating. We did not allow for the possibility that a disoriented 
chiral condensate could appear spontaneously with the nucleation of the bubble. In 
this section we will. For simplicity and definiteness we will use the linear cr-model to 
implement approximate chiral symmetry. The essence of this approach is that the free 
energy difference AF(R) in eq. (2) contains all hadronic degrees of freedom, including 
the pions and the a-meson, but at zero chiral angle. To describe now chiral condensates 
at high temperatures we use the ansatz 

AF(R, 6) = AF(R, 0) + AF chiral (i?, 6) , (15) 

where AF chiral is a perturbative contribution to the free energy difference between hadronic 
and quark-gluon phases for nonzero chiral angle; it is normalized such that AF c hi ra i(-R, 0) = 
0. For low temperatures the system will tend to orient itself towards 9 = 0. On the 
other hand, if the temperature is larger than a certain chiral critical temperature T ch all 
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angles should be approximately equally likely. The aim of this section is to construct 
AF c hi r ai(-R, 9) and to show how to obtain, perturbatively, a distribution in the chiral angle 
over an ensemble of differently-sized disoriented chiral domains. 



3.1 The effective potential for chiral angles 

We will evaluate the finite temperature effective potential representing the degrees of 
freedom of the chiral condensate by using the linear sigma model. 

C = \ (d.af + \ (cVr) 2 - \ (a 2 + vr 2 - c 2 /A) 2 + Ha . (16) 

The parameters of the Lagrangian C are fixed at tree level by imposing the existence of 
spontaneous symmetry breaking, PCAC and the meson masses in vaccuum. This leads 
to the condensates (a) = v and (7Ti) = where 

v{Xv 2 -c 2 )-H = 0, 



2 \ 2 2 

w = Xv -c 

2 _ _2 , o \ „.2 



m 2 = mi + 2\v I . (17) 
Numerical values used here are: /„. = 94.5 MeV, m v = 140 MeV, and m a = 1 GeV. 



The value of the a-mass is chosen in accordance with the analysis of reference |To]. The 
evaluation of the system of equations in (|i~7D yields v = f n , H = 47.6 MeV/fm 2 , A = 54.9 
and c = 686 MeV. 

The finite temperature effective potential for the linear sigma model is discussed in 
many places; for example, in ||16|| . Since eq. (ff^) has a remnant 0(3) symmetry in the 
isovector sector we can simplify this task by letting the condensate point in the third 
direction of the isovector space. To evaluate the thermal masses we expand the fields 
around an arbitrary point as: 

a{x) = v cos 8 + cr'(x) , 

7r 3 (x) = v sin 8 + tt' 3 (x) ■ (18) 

We find 

m l = m 2 = m 3 = ~~ ° 2 ' (19) 

m 2 = 3\v 2 - c 2 . (20) 

In the high temperature limit of the one-loop approximation one keeps terms that are of 
order T 4 and m 2 T 2 only. Then the effective potential is: 

V* (T; v, 8) = ~ [v 4 - (T 2 h - T 2 ) v 2 ] - Hv cos 8 . (21) 
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The chiral critical temperature is defined by T c h = ^2c 2 /X. It is the temperature where 
the curvature term in the effective potential vanishes. When chiral symmetry is exact, 
= 0, and there is a second order phase transition associated with chiral symmetry 
restoration. For our choice of parameters this temperature is 130.9 MeV. It is useful to 
recognize that for a fixed chiral angle 9 the effective potential is minimized when v satisfies 
the cubic equation 

v 3-l(T? h -T 2 )v-^co S 6 = 0. (22) 



2 v cn / A 

We shall label solutions to this equation as v$(T). 

To develop an intuitive understanding of the dynamics of the chiral contribution we 
first plot, in figure 4, the value of the chiral condensate at the global minimum of the 
effective potential as a function of the temperature. The global minimum is always at 
zero chiral angle when chiral symmetry is dynamically broken by nonzero pion mass. 
When chiral symmetry is exact the potential is independent of angle. The dashed curve 
in the figure represents the chiral limit m n = in which we recover a second order chiral 
phase transition; the chiral critical temperature approaches 133.6 MeV when keeping the 
cr-mass and pion decay constant fixed. With the physical value of the pion mass the 
chiral condensate never goes to zero; the minimum of the effective potential is always 
away from the origin and chiral symmetry cannot be restored. Nevertheless the chiral 
condensate v (T) becomes very small at high temperatures, being about 7.9 MeV at 
the confinement temperature of 160 MeV. For high temperatures it decreases like 2 /XT 2 
according to equation (2~2). The minimum of the effective potential itself in ( |2"T| ) decreases 
like 2H 2 / XT 2 . We can now compare the magnitude of the pressure difference of the quark- 
gluon and hadronic phases in eqs. (11)-(12) to the magnitude of the effective potential for 
the chiral angle 9. For more than 2% supercooling the chiral contribution is perturbative. 
There is essentially no nucleation for such small supercooling so the chiral contribution is 
safely in the perturbative region. 

In figure 5 we plot the isotherms of the function vg(T) in the ir^-a plane. At low 
temperatures the isotherms are approximately circles with radii < /„-. The radius shrinks 
and the circle becomes more distorted as the temperature increases. For temperatures 
less than 



T 



saddle 



H\ 2/3 

71-6 (— (23) 



there is a saddle point in the 9 = n direction. This temperature is T sadd i e = 114.9 MeV 
when the pion mass has its physical value. Above this temperature there is no solution 
to the cubic equation at 9 = 0. Rather, the saddle point splits into two and circles back 
around towards the origin with increasing temperature. At a temperature T c h = 130.9 
MeV it reaches the origin. Between T sa ddie and T c h we have a whole range of backward 
angles without a minimum in the radial direction. Above the chiral critical temperature 
we have once again minima in all directions, but they are all centered away from the 
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origin, corresponding to 9 < 7r. Chiral symmetry is approximately restored in the sense 
that vg(T) is very small compared to f n and T. 

The behavior depicted in figure 5 must be dependant on the use of the linear a-model 
to some extent. Nevertheless it is interesting to see that there is a range of temperatures 
where the effective potential does not have minima in the radial direction for parts of 
the backward chiral plane. This would lead to a strong dynamical enhancement of the 
forward direction in the chiral angle distribution. We consider this a typical example of 
the general statement that the inclusion of both dynamical and thermal effects on the 
evolution of distributions tends to smooth out transitions between states and, in this 
sense, suppresses large angle contributions. 

So far we have only discussed the volume term of the chiral angle effective potential. 
Generally one expands the free energy of a finite system in terms of volume, surface, 
curvature and logarithmic contributions. The isotherms of figure 5 are valleys with only 
one minimum located at 9 = 0. Therefore no surface contribution can be defined in 
principle. Only if we include a higher order term in a in the symmetry breaking potential, 
like a quadratic one, can we generate a metastable configuration and therefore define a 
surface free energy. A quadratic term, for example, would generate a barrier at 9 = ±7r/2 
and in this way assure the metastability of the condensate at 9 = it |L6| . 



To construct a surface free energy requires two phases with equal pressures. There 
are, in fact, two phases under consideration, the quark-gluon and the hadronic phases. 
Between these there is a surface free energy (always accepting a first order phase transi- 
tion). However, the 9 dependence of it is not calculable within the linear <j- model, nor 
within any model which does not incorporate the quark and gluon degrees of freedom. 
We take the point of view that the plasma does not care in which direction the chiral field 
points. The 9 dependence of the surface free energy is taken to be zero in this paper. 

3.2 Domain size distributions 

Consider a bubble which was nucleated at a temperature T with a volume V. The 
probability that this bubble was nucleated with a particular chiral angle 9 (constant 
throughout its interior) is 

dV 

— oc exp[-AF chiral (V,T,#)/T] , (24) 

where F^x a \{V, T, 9)/V = V e s(T,vg(T),9) — V e s(T,vo(T),0). If it is assumed that the 
direction of the chiral condensate does not change from the time the bubble was nucleated 
to the end of the phase transition then the joint distribution in chiral angle and domain 
volume at the end of the transition can be computed. From eqs. (13) and (23): 

d 2 n r*t dtt w Nr „ , „ dV 



dVd9 



r l t at t a w 

Jf — I(t)[l-h(t)]8(V-V huh (tM^(V b n h (t,t),T(t),9^ . (25) 
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On the other hand the direction of the chiral condensate may relax according to the 
probability distribution appropriate to the current size and temperature of the bubble. If 
this is the case then the joint distribution at the end of the phase transition would be 

^ = £^/(0[i-MOM(^-W*^))^(W*,*f),r(*f)^) 

= ^ (V, T(tf), 9) £ ^ J(f) [1 - h(t)}5 (v - v hnh (t, u)) 

= %vnw%<y) (26) 

where dn/dV is from eq. (13). It is difficult to know which limit more closely approximates 
reality. The question is: What is the relaxation time for the chiral condensate within a 
bubble? We make an estimate of the critical damping constant which separates the limits 
of weak and strong damping, analogous to a damped harmonic oscillator, in the appendix. 
For the purpose of this paper we compute both limits which will provide bounds on what 
can happen when relaxation effects are included. 

To get a feel for the difference between the extreme limits of underdamping and 
overdamping we plot, in figure 6, the effective potential for the chiral angle divided by 
the temperature versus angle. The Boltzmann/Gibbs factor is obtained by multiplying 
by a volume and exponentiating. The solid curve is for the temperature at the bottom 
of the supercooling curve where the nucleation rate is the highest. The dashed curve is 
for the final temperature at the end of the phase transition. From these curves one can 
already predict that the final distribution in chiral angle will be much broader if there is 
overdamping; that is, if the chiral angle can relax fast enough to follow the temperature 
as it goes up. 

The full double differential distributions in chiral angle and domain size is shown in 
figure 7. Panels (a) and (b) represent parameter sets A and B, respectively, with the chiral 
angle determined at the moment the bubble was nucleated. The "angular distribution" is 
strongly peaked in the forward direction, falling by one to two orders of magnitude as 6 
increases from to tt/2. This is because the dominating nucleation temperature is around 
0.80T C which is comparable to the chiral symmetry breaking scales of v^/V (spontaneous) 
and m n (dynamical). Panels (c) and (d) represent parameter sets A and B, respectively, 
with the chiral angle determined at the end of the phase transition. Since the temperature 
then is greater, about 0.99T C , the angular distribution is much flatter for domain sizes less 
than about 30 V m i n . Thus a large relaxation rate is beneficial to the formation of DCC at 
this point in the heavy ion collision. 

4 Conclusion 

We discussed the dynamical formation of disoriented chiral condensates in ultrarelativis- 
tic heavy ion collisions. Our basic idea is that a first order confinement phase transition 
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governs the dynamical evolution of the system. Once a hot plasma region is formed it 
will cool according to Bjorken hydrodynamics. Bubbles or fluctuations of the low tem- 
perature hadronic phase can form and grow after the system supercools below the critical 
temperature. This process can be described by homogeneous nucleation theory. Every 
bubble or fluctuation can have its own chiral condensate. The orientation of a condensate 
is determined by the Boltzmann/Gibbs weight of the finite temperature effective potential 
of the linear sigma model. We evaluate domain size and chiral angle distributions. 

This model is the first attempt to use the dynamics of a first order phase transition 
to generate disoriented chiral condensates. Many of the approximations we made in 
this model can be relaxed in more elaborate calculations. For example, the influence of 
the chiral contribution to the free energy is treated perturbatively to lowest order. The 
dependance on rapidity, the eccentricity of the domains, is neglected. The nucleation 
description we apply is insufficient both at the very early stages of the supercooling 
process, where very large bubbles containing several different domains might form, as 
well as in the final stages of the phase transition, where surface and/or topological effects 
between overlapping domains might become important. An improved description of the 
thermodynamic potential for the chiral degrees of freedom as well as an investigation of 
the importance of dissipation in the evolution of a chiral domain embedded in a hadronic 
heat bath is desirable. 

It is currently accepted that both the confining character of QCD as well as the 
properties imposed by an approximate chiral symmetry are crucial in determining the 
properties of the QCD phase transition. In the model presented we intended to incorporate 
both properties in a sensible description of the phase transition dynamics and the resulting 
collective observables. We are of the opinion that modifications and improvements of our 
model as well as different approaches in the same class of models will be crucial in not 
only addressing issues like the formation of DCCs but also in addressing all issues related 
to the late stages of the quark-hadron phase transition. 
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Appendix 



In this appendix we study small amplitude oscillations of the chiral angle about the ab- 
solute minimum of the effective potential to a get an estimate of the time scales involved 
relative to the time evolution of the temperature. Replacing the potential in the La- 



grangian (16) with the effective potential in (21) results in the equation of motion for 
9. 

v 2 d 2 9-Hvsm9 = 0. (27) 

This is the same equation one would have obtained by neglecting the thermal one-loop 
corrections. However, it does not include thermal effects on the derivative terms. The 
thermal effects are incorporated via the radial field v which will be constrained to satisfy 
dV e g/dv = in equation (|2^) . Hence we replace v in eq. fl27|) with vo(T) thus eliminating 
its thermal fluctuations. 

We would like to solve eq. (^) deep inside an hadronic bubble where we can neglect 
any spatial gradients. The equation of motion conserves energy and therefore cannot 
describe the decay of 9 into the absolute miniumum, which must happen if the temperature 
decreases to zero. To obtain a simple order of magnitude estimate of the time scale in the 
decay process we introduce a dissipative term of the form r]v 2 d9/dt. In this ansatz 7] has 
the dimension of energy and corresponds to a Langevin type of friction term. For small 
oscillations around zero, relevant at low temperatures: 

v 2 (T) 6 + VVq(T) 9 + Hv (T) 9 = 0. (28) 

The oscillation frequencies uj± are 




UJ± = ^ ± ^~ r]2 + ^{T)- (29) 



The oscillator is critically damped for the value 



Vc = 2 J-^r 2m* , T < T ch . (30) 
V M T ) 

Since this is a typical hadronic scale it is not clear whether the chiral angle will be under- 
or over-damped. A stronger damping than this allows the system to relax very rapidly 
on time scales of less than r/^ 1 . 

For temperatures larger than T c ^ we see from figure 5 that the isotherms are approx- 
imately circles in the forward half-plane with radius l/2t> (T) centered at (l/2t> (T), 0). 
Expanding eq. (|2"7| ) around the center of this circle we retain eqs. (^) and ( PUD but with 
Vq and scaled by 1/2. Using the results from section 2 about the high temperature 
behavior we find that 

r/ c = VAT, T>T ch . (31) 
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The friction constant would have to increase linearly with the temperature at high tem- 
peratures to remain over-damped. This result can be understood as a consequence of 
approximate chiral restoration at high temperatures forcing vq(T) to become small. Ac- 
tually, our small angle approximation becomes ill-defined then. 
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Set A 


Set B 




4.6 


7.5 


g q 


37 


47.5 


Bi (MeV) 


220 


232 


t sc (fm/c) 


6.7 


7.2 


T sc /T c 


0.81 


0.79 


t{ (fm/c) 


30.9 


25.4 


T f /T c 


0.991 


0.987 


£ min (fm/c) 


8.7 


9.7 


K„in (fm 3 ) 


21.8 


10.1 


a 


0.67 


0.41 



Table 1. Summary of observables and scales for the two different parameter sets. 
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Figure Captions 



Figure 1. The time evolution of the temperature. 

Figure 2. The bubble volume at the end of the phase transition as a function of the time 
it was nucleated. 

Figure 3. The domain size distribution at the end of the phase transition. 
Figure 4. The thermal average value of the <r— field. 

Figure 5. Polar plot of the radial minimum of the finite temperature effective potential 
in the n 3 -a plane. 

Figure 6. The cost in free energy per unit volume divided by the temperature as a 
function of the chiral angle. The solid curve is evaluated at the moment the bubble is 
nucleated while the dashed curve is evaluated at the end of the phase transition. 

Figure 7. The double distribution V^ in d 2 n/dV d9 in chiral domain size and direction at 
the end of the phase transition. In panels (a) and (b) the direction is fixed at the moment 
of bubble nucleation and in panels (c) and (d) the direction is fixed by the temperature 
at the end of the phase transition. 
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